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Water and Methanol Maser Activities in the NGC 2024 FIR 6 Region 
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ABSTRACT 

The NGC 2024 FIR 6 region was observed in the water maser line at 22 GHz and the methanol class I maser 
lines at 44, 95, and 133 GHz. The water maser spectra displayed several velocity components and month-scale time 
variabilities. Most of the velocity components may be associated with FIR 6n, while one component was associated with 
FIR 4. A typical life time of the water-maser velocity-components is about 8 months. The components showed velocity 
fluctuations with a typical drift rate of about 0.01 km s _1 day -1 . The methanol class I masers were detected toward FIR 
6. The methanol emission is confined within a narrow range around the systemic velocity of the FIR 6 cloud core. The 
methanol masers suggest the existence of shocks driven by either the expanding H II region of FIR 6c or the outflow of 
FIR 6n. 

Subject headings: ISM: individual (NGC 2024 FIR 6) — ISM: structure — Masers — stars: formation 



1. INTRODUCTION 

The NGC 2024 region is a site of star formation at a dis- 
tance of 415 pc from the Sun (Anthony-Twarog 1982). FIR 6 
is one of the dense cores in the NGC 2024 molecular ridge 
(FIR 1-7), first identified by observations in the dust continuum 
(Mezger et al. 1988), and was known to be associated with an 
H 2 maser source and a molecular outflow (Genzel & Downes 
1977; Richer 1990). Studies with the millimeter continuum and 
molecular lines suggested that FIR 6 may contain a protostar 
(Schulz et al. 1991; Chandler & Carlstrom 1996; Wiesemeyer 
et al. 1997; Visser et al. 1998; Mangum et al. 1999). 

High-resolution imaging observations revealed a binary sys- 
tem (Lai et al. 2002). The projected separation between the bi- 
nary members, FIR 6c and 6n, is 3'.'9 or 1600 AU, and both of 
them are associated with H2O maser sources (Choi et al. 2012). 
The primary (in the millimeter continuum), FIR 6c, has a flat 
continuum spectrum around 1 mm and suggested to be a hy- 
percompact H II region probably powered by a Bl star (Choi et 
al. 2012). The H2O maser line of FIR 6c detected in 1996 was 
blueshifted by ^20 km s _1 with respect to the ambient cloud 
velocity, suggesting the existence of an outflow activity (Choi 
et al. 2012). The secondary, FIR 6n, exhibits relatively stronger 
H9O maser emission and drives a bipolar outflow (Furuya et al. 
2003; Alves et al. 2011; Choi et al. 2012). FIR 6n may be a 
low-mass young stellar object (YSO), probably an active proto- 
star (Choi et al. 2012). As the observed phenomena displayed 
by FIR 6 are a mixture of high- and low-mass star formation 
activities, further investigations are needed to understand this 
interesting region. 

In this paper, we present the results of our observations of the 
NGC 2024 FIR 6 region in the H 2 0, CH 3 OH, and 13 CS lines 
with the Korean Very Long Baseline Interferometry Network 
(KVN) antennas. We describe our spectral line observations in 
Section 2. In Section 3, we report the results. In Section 4, we 
discuss the star-forming activities in the FIR 6 and 4 regions. A 
summary is given in Section 5. 

2. OBSERVATIONS 

The NGC 2024 FIR 6 region was observed using the KVN 2 1 
m antennas in the single-dish telescope mode during the 201 1— 
2012 observing season. The observations were carried out with 
the KVN Yonsei telescope at Seoul, Korea in 2011 November 



and 2012 January, and the KVN Ulsan telescope at Ulsan, Korea 
in 2012 April-May. 

The data were obtained using the high-electron-mobility- 
transistor receivers in the 22, 44, and 86 GHz bands and the 
superconductor-insulator-superconductor receivers in the 129 
GHz band. Digital filters and spectrometers were used as the 
back end (Lee et al. 2011). The spectrometers were used 
in the mode of 4 intermediate-frequency data streams. Each 
intermediate-frequency stream was configured to have a band 
width of 64 MHz and 4096 channels, which gives a spectral 
channel width of 15.625 kHz. The antenna temperature was cal- 
ibrated by the standard chopper-wheel method, which automat- 
ically corrected for the effects of atmospheric attenuation. 

The main target source was NGC 2024 FIR 6 at 02000 
= 05 h 41 m 45!14, <5 200 o = -0r56'04'/4, which corresponds to 
the 6.9 mm continuum source FIR 6c (Choi et al. 2012). The 
data were taken by position switching with the off position of 
(05 h 43 m 05!19, -01°56'04'.'4), which is a nearby empty area in 
the CS J = 2 -» 1 map of Lada et al. (1991). 

Telescope pointing was checked by observing Orion IRc2 
(Baudry et al. 1995) in the SiO u=lJ=l->0 maser line. The 
pointing observations were performed about once in an hour. 
The telescope pointing was good within ~4". 

2.1. Main Target Lines 

The main target lines were the H 2 616 — > 593 (22 GHz) line 
and the CH 3 OH 7 -> 61 A+, 8 -> 7 a A+, and 6_ x -> 5 E 
(44, 95, and 133 GHz, respectively) lines. These CH 3 OH lines 
are class I maser lines. Table 1 lists the parameters of these lines. 
The spectrometer setting gives a velocity channel width of 0.21 
km s- 1 for the H 2 line and 0.11 km s" 1 for the CH 3 OH 44 
GHz line, and we present the data of these lines without fur- 
ther resampling. For the CH3OH 95 GHz line, the spectra were 
smoothed using a Hanning window over three channels, which 
gives a velocity channel width of 0.098 km s _1 . For the CH3OH 
133 GHz line, the Hanning smoothing was applied twice, which 
gives a velocity channel width of 0. 14 km s _1 . 

In the 2012 May run, a map was made in the H 2 line, cov- 
ering a 4.'8 x 4.'8 region around the FIR 1-7 cores. The map- 
ping grid size was 12" near FIR 3-7 and 48" in the periphery 
of the map. Some extra integrations were made toward FIR 4 at 
(05 h 41 m 44!14, -01°54'46"0). 

The data were calibrated using the standard efficiencies of 
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Table 1 

Parameters of Detected Maser Lines 







Frequency a 




Line 




(GHz) 


Observing Run 


H 2 6ie -s 


' 523 


22.2350771 


2011 Nov 6, 2012 Jan 27, 2012 Apr 12-14, 2012 May 11-12 


CH3OH 7 


->6i A+ 


44.069476 


2012 Apr 12-15, 2012 May 1 1-13 


CH3OH 80 


->■ 7i A+ 


95.169516 


2012 Jan 27, 2012 Apr 12-15 


CH3OH 6_ 


1 -+ 5 E 


132.890800 


2012 Apr 12-15 



a Rest frequency from the line catalogue of the National Institute of Standards and Technology (NIST) 
recommended rest frequencies for observed interstellar molecular microwave transitions by F. J. Lovas 
(http://www.nist.gov/pml/data/micro/index.cfm). 



Table 2 
Telescope Parameters 







Frequency 3 


Beam b 






f 


Telescope 


Observing Run 


(GHz) 


(arcsec) 




?7mb C 


(Jy K- 1 ) 


KVN Yonsei 


2011 Nov, 2012 Jan 


22 


118 


0.65 


0.46 


15.3 






86 


28 


0.48 


0.36 


20.7 


KVN Ulsan 


2012 Apr, 2012 May 


22 


119 


0.62 


0.44 


16.0 






44 


62 


0.62 


0.47 


16.0 






86 


29 


0.49 


0.39 


20.3 






129 


22 


0.32 


0.29 


38.8 



a Frequency where the beam size and efficiencies were measured. 
b Full width at half maximum (FWHM) of the main beam. 

c Aperture and main-beam efficiencies. In each receiver band, t)a is nearly constant, and i]mb 
is inversely proportional to the square of frequency. 

d Scaling factor for converting the KVN raw data to spectra in the flux density scale for 
the maser lines in Table 1, which includes the telescope sensitivity, quantization correction 
factor (1.25), and sideband separation efficiency (0.8 for the KVN Ulsan 129 GHz band and 1.0 
otherwise). 



KVN for the 2011-2012 observing season (Table 2; also see 
Lee et al. 2011 and http://kvn-web.kasi.re.kr/). We will present 
the maser data in the flux density scale. For each spectrum, a 
first-order baseline was removed. The baseline was determined 
from the velocity intervals of Vlsr = (-10, 0) and (20, 30) km 
s _1 for the H2O line and from (1, 6) and (16, 21) km s _1 for 
the CH3OH lines. 

2.2. Other Methanol Lines 

In addition to the class I maser lines described above, FIR 6 
was observed in most of the known CH3OH lines in the tuning 
ranges of the KVN receivers. These observations were carried 
out in the 2012 April-May observing runs. Table 3 lists the line 
parameters. None of these lines was detected. 

2.3. The 13 CS J = 2 1 Line 

FIR 6 was also observed in the 13 CS J = 2 — > 1 line 
(92.494270 GHz) in the 2012 April observing run. The spec- 
tra were smoothed using a Hanning window over three chan- 
nels, which gives a velocity channel width of 0.10 km s _1 . For 
each spectrum, a first-order baseline was removed. The spec- 
tral baseline was determined from the velocity intervals of (1, 
6) and (16, 21) km s -1 . We will present the 13 CS data in the 
main-beam temperature (T m b) scale. 



3. RESULTS 

3.1. The Water Maser 

The H2O maser line was detected in all the four observing 
runs. The spectra showed several velocity components, and each 
of them displayed significant month-scale time-variability. Fig- 
ure 1 shows the spectra from each observing run. 

Inspecting the spectra closely, five velocity components can 
be identified: Vlsr ~ 3, 5, 6, 8, and 13 km s _1 . The 3 km 
s _1 component was present in the first observing run, became 
stronger in 2012 January, weakened in April, and then disap- 
peared. The 5 km s _1 component was present in 2012 April 
(as a shoulder of the 6 km s _1 component), and could be seen 
clearly in May. The 6 km s _1 component "flared" in 2012 April, 
and weakened in May. The 8 km s _1 component appeared in 
2012 April, and became stronger in May. The 13 km s _1 com- 
ponent appeared in 2012 January, became stronger in April, and 
weakened in May. The map data showed that most of the ve- 
locity components are associated with FIR 6 while the 13 km 
s _1 component is associated with FIR 4 (Section 3.1.1). Table 
4 lists the line-profile parameters (centroid velocity, peak inten- 
sity, line FWHM, and integrated intensity) derived by Gaussian 
fits to each velocity component. 
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Table 3 

Parameters of Undetected Methanol Lines 







Frequency a 


r — 

rms 


rms c 


Line 




(GHz) 


(K) 


(Jy) 


CH3OH 12 2 


->■ lli A~ v t = 1 


21.550342 


0.16 




CH3OH 9 2 - 


-»■ 10i A+ 


23.121024 


0.14 


0.7 


CH3OH 6-2 


-»■ 7-i E 


85.568074 


0.27 


1.7 


CH3OH 7 2 - 


-> 63 A" 


86.615602 


0.24 


1.5 


CH3OH 7 2 - 


-> 63 


86.902947 


0.21 


1.3 


CH3OH 15 3 


->■ 144 A+ 


88.594809 


0.23 




CH3OH 15 3 


-> 144 A - 


88.939993 


0.21 




CH3OH 8-4 


-> 9-3 £ 


89.505778 


0.31 




CH3OH 20- 


3 -> 19-2 Ev t = l 


90.81239 


0.29 




CH3OH 83 - 


-y 9 2 £ 


94.541806 


0.31 


1.6 


CH3OH 60 - 


-> 5i E 


124.569976 


0.45 




CH3OH 12i 


-> 11 2 A" 


129.433406 


0.41 




CH3OH 62 -> 7i A - 


132.621859 


0.43 


3.0 


CH3OH 5-2 


-> 6-1 £ 


133.605385 


0.87 


5.9 


CH3OH 12_ 


3 -> 13-2 £ 


134.231013 


0.37 




CH3OH 7-4 


-»• 8-3 E 


137.902997 


0.41 




CH3OH 23- 


2 -»■ 23i E 


140.03314 


0.57 




CH3OH 18c, -)■ I8-1 £ 


140.151188 


0.38 




CH3OH 0o - 


■> li Ev t = 1 


141.441249 


0.52 





a Rest frequency from the NIST catalogue. 

b Noise rms level in the T m b scale, measured from spectra with a 
channel width of 15.625 kHz (0.2 km s" 1 ) for the 22 GHz band, 
31.25 kHz (~0.11 km s" 1 ) for the 86 GHz band, and 62.5 kHz 

(~0.14 km s _1 ) for the 129 GHz band. 

c Noise rms level in the flux density scale for class II maser lines 
(Menten 1991; Cragg et al. 1992; Miiller et al. 2004). 



Table 4 

Parameters of the Water Maser Spectra 







Sp 


AV a 


J" Sdv 


Run 


(km s^ 1 ) 


(Jy) 


(kms^ 1 ) 


(Jykms- 1 ) 


2011 Nov 


3.34 


5.1 


0.72 


3.9 ± 0.7 


2012 Jan 


2.74 


8.0 


0.84 


7.2 ± 0.3 




13.41 


11.1 


1.03 


12.1 ±0.3 


2012 Apr 


2.98 


4.1 


1.17 


5.1 ±0.5 




5.38 


21.1 


1.85 


41.5 ±5.4 




6.38 


63.1 


1.32 


88.8 ±5.4 




8.09 


16.9 


0.64 


11.4 ±0.3 




12.86 


32.4 


0.71 


24.4 ± 0.3 


2012 May b 


4.91 


42.1 


0.77 


34.7 ± 1.3 




6.19 


11.5 


1.06 


12.9 ± 1.6 




7.57 


43.0 


0.88 


40.4 ± 1.3 




13.17 


49.4 


0.89 


46.6 ± 0.9 



a Uncertainties of the centroid velocity and line FWHM are 
much smaller than the channel width (0.21 km s _1 ). 

b The parameters of the 5, 6, and 8 km s _1 components are 
from the spectrum toward FIR 6 (Figure 1(d)), and those of the 
13 km s _1 component are from the spectrum toward FIR 4 (Fig- 
ure 1(e)). 
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Fig. 1. — Spectra of the H2O 22 GHz maser line, (a-d) Spectra toward 
FIR 6. (e) Spectrum toward FIR 4. The observation dates are labeled. 
Vertical dotted line: velocity of the ambient dense gas (Vlsr. = 10.8 
km s _1 ; Section 3.3). 



3.1.1. Source Positions 

Examinations of the map data showed that the 5, 6, and 8 
km s _1 components come from sources near FIR 6, and they 
(if different) cannot be distinguished with the KVN beam. By 
contrast, the 13 km s _1 component clearly comes from a differ- 
ent source near FIR 4. (The 3 km s _1 component had already 
disappeared, and its relation to these sources is not clear.) Fig- 
ure 2 shows the H2O line maps. The intensity distribution in 
each panel is consistent with what is expected from a point-like 
source (convolved with the beam). To find the source position, 
each map was fitted with a Gaussian intensity profile having the 
same FWHM as the beam. 

The best-fit source position of the 5/6/8 km s _1 components 
is (-2" ± 7", 7" ± 6") with respect to FIR 6c. Among the 
two known H2O maser sources in this region (Choi et al. 2012; 
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Fig. 2. — Maps of the H2O 22 GHz maser emission in 2012 May. The 
contour levels are 25, 50, and 75% of the maximum value in each panel, 
(a) Map of the 5/6/8 km s" 1 components, (b) Map of the 13 km s _1 
component. The velocity interval for the integration of intensity is la- 
beled in each panel. Plus signs: H2O maser source positions derived 
by Gaussian fits. The size of plus sign corresponds to the position un- 
certainty. Diamond symbols: millimeter continuum sources FIR 1-7 
(Chandler & Carlstrom 1996; Visser et al. 1998; Choi et al. 2012). 



Furuya et al. 2003), FIR 6n is located within the position un- 
certainty boundary, and FIR 6c is just outside the boundaryQ 
We presume that the source of the 5/6/8 km s _1 components is 
closely associated with FIR 6n, while FIR 6c cannot be com- 
pletely ruled out. 

3 As the description of the H2O maser in Furuya et al. (2003) is prone to con- 
fusion, we reanalyzed the data set of their observations with Very Large Array 
(project AF 354). We found that the source coordinate in their Table 3 is cor- 
rect, and it corresponds to FIR 6n. No H2O maser emission was detected in the 
region around FIR 5. The source described as "7000 AU southeast of the close 
binary FIR 5" in their Section 4.3.4 should be interpreted as FIR 6n. 
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The best-fit source position of the 13 km s _1 component is 
(-1" ± 9", 4" ± 8") with respect to FIR 4. There is no previous 
report of an H2O maser detection in this region. Since FIR 4 is 
the only known YSO within the uncertainty boundary, we pre- 
sume that the source of the 13 km s _1 component is associated 
with FIR 4. 

Note that the first detection of the H2O maser in the Orion 
B region was reported by Johnston et al. (1973). However, their 
position uncertainty was large (1'), and FIR 3-6 are in the un- 
certainty boundary. Therefore, it is difficult to identify the object 
responsible for the H 2 emission detected by them. 

3.1.2. Variabilities 

Month-scale variabilities of the H2O maser flux density are 
shown in Figure 3(a). Since four of the velocity components 
showed flux maxima in the middle observing runs, a rough es- 
timate of timescale can be derived. The FWHM of each light 
curve ranges from 60 to 180 days. If the occurrence of each 
velocity component is a one-off event (that is, shows a peak 
flux and disappears without a recurrence), a typical life time 
(2FWHM) would be -8 months. 

Figure 3(b) shows variabilities of the centroid velocity. Inter- 
estingly, in all cases, the velocity drifted to a lower value (that 
is, blueshifted) in the first time interval after the appearance. It 
is unlikely that this phenomenon is owing to an instrumental 
artifact because the amount and sense of the drifts are all dif- 
ferent even in the same time interval. (For example, in the last 
time interval of Figure 3(b), the 13 km s _1 component drifted 
to a higher velocity while the others drifted to lower velocities 
at various rates.) Those detected three times (the 3 and 13 km 
s _1 components) changed the sense of velocity drift. Therefore, 
the centroid velocity is fluctuating rather than drifting steadily. 
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Fig. 3. — Variabilities of the H2O 22 GHz maser for each velocity com- 
ponent. The horizontal axis shows the time since the first observing 
run. (a) Light curve. For the 3, 5, 6, and 8 km s _1 components, the in- 
tegrated flux densities are from the spectra toward FIR 6. For the 13 km 
s _1 component, the integrated flux densities of the 2012 January /April 
runs are from the spectra toward FIR 6, corrected for the primary beam 
attenuation, and that of the 2012 May run is from the spectrum to- 
ward FIR 4. (b) Time evolution of centroid velocity with respect to the 
velocity at the first detection. The velocities of the identified velocity 
components are labeled. 



A typical drift rate is ^0.01 km s _1 day -1 . The velocity drift is 
usually interpreted as the acceleration of emission source (e.g., 
Tofani et al. 1995), but the fluctuations of the velocity seen in 
the FIR 4/6 masers suggest that the cause of velocity variability 
is more complicated. 

Day-scale variabilities were studied in the 2012 April ob- 
serving run. Figure 4 shows the change of spectra in two days. A 
difference of line profile can be seen in the 6 km s _1 component 
only. This component became stronger by 10% and redshifted 
by ^0.09 km s _1 . (For comparison, the intensity variation of the 
13 km s _1 component was smaller than 2%.) The information 
we can extract from the two spectra is rather limited. Neverthe- 
less, it is clear that the sign of velocity drift in these two days 
is the opposite of what was seen in a month as describe in the 
paragraph above, which reinforces the finding that the velocity 
of emission source fluctuates. 

3.2. The Methanol Masers 

All the three CH3OH class I maser lines in the tuning ranges 
of the KVN receivers were detected toward FIR 6. By contrast, 
other CH3OH lines were undetected (Table 3), which suggests 
that the detected lines are amplified emission. Each of the 44 
and 95 GHz lines was observed in two observing runs, and there 
was no detectable time variability. 

Figure 5 shows the spectra for each line, and Table 5 lists 
the line-profile parameters. All the detected CH3OH emission 
is confined in a velocity interval around the systemic velocity, 
within ^2 km s _1 . For the 44 GHz line, the residual spectra 
of single-component and double-component Gaussian fits show 
a peak at ~10 km s _1 higher than the noise level. Therefore, 
the 44 GHz line seems to show three velocity components. The 
signal-to-noise ratios of the 95 and 133 GHz lines are too small 
to tell if there are multiple velocity components. 

Considering the non-detection of the other CH3OH lines and 
the small width of each velocity component of the 44 GHz line, 
we favor the interpretation that the detected lines are maser 
emission. However, it is possible that the detected emission can 
be a mixture of thermal and maser components (e.g., Kalen- 
skii et al. 2010). Additional observations such as interferometric 
imaging are needed to confirm the nature of this emission. 
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Fig. 4. — Daily spectra of the H2O 22 GHz maser line toward FIR 6 in 
the 2012 April run. (a) Spectra on April 12 (solid line) and 14 (dotted 
line), (b) Difference between the two days. 




5 



Table 5 

Parameters of the Methanol Maser Spectra 





V 


Q 

Op 


AT/ 

Av 


J SdV 


Line 


n™, n -l\ 

(Km s ) 


(Jy) 


(km s ) 


(Jy Km s J 


CH3OH 44 GHz a 


11.12 ±0.08 


1.4 


2.2 ± 0.2 


3.2 ±0.2 


CH3OH 44 GHz b 


10.32 ± 0.08 


1 ^ 


U.o It u.z 


1.1 It KJ.J 




11.32 ±0.08 




0.7 ± 0.2 


1.2 ± 0.4 




12.21 ±0.14 


0.9 


0.7 ± 0.3 


0.7 ± 0.3 


CH3OH 95 GHz a 


10.99 ±0.13 


2.5 


1.8 ±0.3 


4.7 ± 0.7 


CH3OH 133 GHz a 


10.41 ±0.19 


7.7 


1.7 ± 0.4 


13.7 ±2.7 


a Single-component Gaussian fit. 








b Triple-component Gaussian fit. 
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Fig. 5. — Spectra of the CH3OH class I maser lines toward FIR 6. (a) 
CH3OH 7 61 A+ line, (b) CH 3 OH 8 -> 7i A+ line, (c) CH3OH 
6-1 -r 5o line. For each line, all the data from the 2011-2012 
season were averaged. 



22 (Class I YSO; Comeron et al. 1996) is located within the 
beam. 

3.3. The 13 CS J = 2^1 Line 

The 13 CS line was observed to measure the systemic velocity 
of the dense cloud core containing FIR 6. Figure 6 shows the 
spectrum. The best-fit Gaussian profile gives Vq = 10.77 ± 0.05 
km s"\ peak T mb = 1.95 ±0.03 K, AV = 1.66 ±0.11 km s~\ 
and jT mh dV = 3.44 ± 0.19 K km sM 

The systemic velocity of the FIR 6 dense core has been 
known somewhat poorly because the values in the literature 
were measured using lines with relatively higher optical depths. 
Schulz et al. (1991) derived 11.3 km s" 1 using the CS/C 34 S 
J = 5 — > 4 and 7^6 lines showing widths of 2.5-2.9 km 
s _1 . Mangum et al. (1999) derived 10.6-10.9 km s _1 using the 
H2CO lines showing widths of 2.3-2.5 km s _1 . For the deriva- 
tion of systemic velocity, 13 CS is a better tracer than C 34 S be- 
cause of its smaller abundance (Shirley et al. 2003). Indeed, the 
13 CS J = 2 — > 1 line is much narrower than the lines listed 
above. Therefore, the systemic velocity of 10.8 km s^ 1 derived 
from the 13 CS line may be more reliable. 

4. DISCUSSION 
4.1. FIR 6 

The detection of the H2O and CH3OH masers suggests that 
the NGC 2024 FIR 6 region is actively forming stars, which 
produces shocks in the warm and dense molecular gas around 
the YSOs. The two maser species show the well-known con- 
trasting characteristics (Menten 1991; Elitzur 1992). The H 2 



There is no map of the CH 3 OH lines with a sensitivity good 
enough to constrain the source position, and the YSO(s) respon- 
sible for the maser emission can only be deduced from the beam 
size. Since these lines have different beam sizes (Table 2), it 
is possible that some emission structure in the beam of a low- 
frequency line can be outside the beam of a high-frequency line. 
It is difficult to specify the source of the 44 GHz maser line be- 
cause there are several embedded objects (FIR 5w/e, 6c/n, and 
7) and other YSOs within the beam (Barnes et al. 1989; Skinner 
et al. 2003; Rodriguez et al. 2003; also see Figure 1(b) of Choi 
et al. 2012). For the 95 and 133 GHz maser lines, the source is 
most likely associated with FIR 6c/n because the angular sepa- 
rations of FIR 5 and 7 from FIR 6c ( -25 " and 2 1 " , respectively) 
are larger than or comparable to the beam sizes. Note that IRS 
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Fig. 6.— Spectrum of the "CS J = 2 -> 1 line toward FIR 6. Vertical 
dotted line: centroid velocity of the best-fit Gaussian profile (Vlsr = 

10.8 km sM. 
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maser emission can be seen in a relatively wide range of veloc- 
ities while the CH3OH maser emission is confined in a narrow 
range around the systemic velocity of the cloud. The H2O maser 
also displays strong variabilities in both amplitude and velocity, 
while the CH3OH masers show no detectable variability. 

The most likely exciting source of the H2O maser reported in 
this paper is FIR 6n. The maser probably comes from the base 
of the bipolar outflow driven by FIR 6n (Richer 1990; Choi et 
al. 2012). During the observation period covered in this study, 
all the H 2 emission was blueshifted with respect to the sys- 
temic velocity of the FIR 6 core, and the blueshifted outflow 
may have been responsible for the observed H 2 maser emis- 
sion. However, the emission was nearly at the systemic velocity 
in 1995/1996 (Choi et al. 2012), and both blueshifted and red- 
shifted velocity components were detected in 1999 (Fumy a et 
al. 2003). 

The variability of the FIR 6 H 2 maser is typical of low- 
mass YSO. The typical life time of each velocity component 
is —8 months. This timescale is generally consistent with the 
values reported in surveys toward low-mass YSOs, which range 
from —2 weeks to —8 months (Wilking et al. 1994; Furuya et 
al. 2003). The velocity drift rate of the FIR 6 H 2 maser is 
similar to that reported by Tofani et al. (1995). Note that the 
time interval between observing runs in this study is not short 
enough to trace variabilities at timescales shorter than about a 
month. 

The interpretation of the CH3OH class I masers is unclear. 
Soon after the discovery of the CH 3 OH masers, it was found 
that class I masers trace outflows driven by high-mass YSOs 
and are located away from the responsible YSOs up to —2 pc 
(Plambeck & Menten 1990; Kurtz et al. 2004). Recently, how- 
ever, it was found that some class I masers may be associated 
with shocks driven into molecular clouds by expanding H II re- 
gions (Voronkov et al. 2010; Chen et al. 201 1) and that several 
low-mass star forming regions emit class I masers (Kalenskii et 
al. 2010). If the FIR 6 CH3OH masers trace an expanding H II 
region, the FIR 6c hypercompact H II region (candidate) may be 
responsible. If they trace an outflow, the FIR 6n bipolar outflow 
may be responsible for the masers. If the maser emission comes 
from the FIR 6n outflow, it would be one of the rare examples 
of CH3OH maser associated with low-mass YSOs. This issue 
can be settled by future observations using interferometers. 

4.2. FIR 4 

The detection of the H 2 maser toward FIR 4 illustrates 
the importance of mapping in single-dish observations of H 2 
masers. While this paper gives the first confirmation of the H 2 
maser associated with FIR 4, it is not necessarily the first de- 
tection. For example, Furuya et al. (2003) reported single-dish 
observations toward FIR 5/6 (see their Figure 13(a)). Since FIR 
4-6 are within their beam, it is likely that the detected emis- 
sion in their spectra were mostly from FIR 4/6. (Note that in- 
terferometric observations never showed an H 2 maser associ- 
ated with FIR 5.) Our recent survey also confirmed that multiple 
H 2 maser sources in a single-dish beam is a common occur- 
rence in cluster-forming regions such as the Orion cloud (M. 
Kang et al. 2012, in preparation). The exciting source of H 2 
maser emission in a given detection is not necessarily the pre- 
viously known source, the youngest YSO in the beam, or the 
most powerful outflow source in the beam. This problem can 
affect some statistical quantities derived from single-dish sur- 
veys, such as detection rates, and can be worse in massive star 
forming regions far away from the Sun. 

The H 2 maser suggests that there are shocks produced by 



the star formation activities of FIR 4. The nature of FIR 4, how- 
ever, has been controversial. Early studies indicated the exis- 
tence of dense gas in the NGC 2024 region (Evans et al. 1987). 
Mezger et al. (1988, 1992) suggested that the FIR 1-7 cores are 
isothermal protostars without luminous stellar objects, based 
on their dust continuum observations. Studies with molecular 
(thermal) lines suggested that they contain more evolved YSOs 
(Schulz et al. 1991; Mangum et al. 1999). The discovery of a 
near-IR reflection nebula and a CO outflow revealed that FIR 4 
contains a low-mass YSO driving a bipolar outflow (Moore & 
Chandler 1989; Moore & Yamashita 1995; Chandler & Carl- 
strom 1996). A very different view was proposed by Minier 
et al. (2003). They detected the CH3OH class II maser line at 
6.7 GHz in emission and argued that FIR 4 contains an inter- 
mediate or high-mass YSO. Future investigations are needed to 
tell whether low-mass protostars can produce CH3OH class II 
maser emission or FIR 4 is a site of high-mass star formation. 
The detection of the H 2 maser does not necessarily favor one 
or the other possibility. 

It is interesting to note that there are position offsets of 
—2" among the 3 mm continuum, CH3OH class II maser, and 
hard X-ray positions (Chandler & Carlstrom 1996; Minier et al. 
2003; Skinner et al. 2003). Future determinations of the accu- 
rate position of the H 2 maser source will be helpful in under- 
standing the nature of FIR 4. 

5. SUMMARY 

The NGC 2024 FIR 6 region was observed using the KVN in 
the H 2 22 GHz maser line, the CH 3 OH 44, 95, and 133 GHz 
class I maser lines, other CH3OH lines in the tuning ranges, and 
the 13 CS J = 2 — > 1 line in the single-dish telescope mode. In 
addition, the FIR 4 region was observed in the H 2 line. The 
main results are summarized as follows. 

1. Several velocity components were identified in the H 2 
spectra. The H 2 maps revealed that most of them seem to be 
associated with FIR 6n while one is associated with FIR 4. None 
is associated with FIR 5. 

2. Month-scale variabilities of the H 2 maser were seen 
clearly. A typical life time of each velocity component is —8 
months. The centroid velocity fluctuates with a drift rate of 
-0.01 km s" 1 day" 1 . 

3. The CH3OH class I maser lines were detected toward FIR 
6. They were interpreted as maser emission. The maser emis- 
sion is confined within —2 km s -1 around the systemic veloc- 
ity of the cloud. It is not clear whether the exciting source is 
FIR 6c or 6n. The CH 3 OH masers did not show a detectable 
time-variability. High-resolution imaging is needed to under- 
stand whether the CH3OH maser emitting region is shocked by 
an outflow or an expanding H II region. 

4. The other CH 3 OH lines, thermal or class II maser, were 
not detected. 

5. The 13 CS spectrum showed that the systemic velocity of 
the FIR 6 dense core is —10.8 km s -1 . 
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